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1 Introduction 

Recent work by Janik and collaborators [1, 2, 3] (see also [4]) opened a possibility 
to study non-stationary/non-equilibrium processed in gauge theories using their dual 
string theory formulation [5, 6]. In fact, one might try to use gauge theory/string 
theory correspondence to define non-equilibrium Quantum Field Theory dynamics 1 . 
For example: in principle, one can use AdS/CFT correspondence to extend Muller- 
Israel-Stewart theory of dissipative relativistic hydrodynamics [8, 9] to arbitrary order 
in deviation from the equilibrium. 

Even though the framework proposed in [1] can not be consistently implemented 
within a supergravity approximation of the string theory dual to M = 4 supersymmet- 
ric Yang-Mills (SYM) plasma in the Bjorken flow [10] (see [11]), if correctly reproduces 
the equilibrium and near-equilibrium properties of the M = 4 SYM plasma obtained 
from the equilibrium correlation functions, i.e., the equilibrium equation of state [12], 
its shear viscosity at infinite 't Hooft coupling [13], and its relaxation time 2 [17]. Fur- 
thermore, this agreement appears to be unaffected by the subtleties associated with 
the breakdown of the supergravity approximation. Given substantial potential of the 
Janik's framework 3 , we believe it is important to further explore its regime of validity. 

As we already mentioned, the shear viscosity of the M = 4 plasma at infinite 
coupling is a robust prediction within Janik's framework [2]. In this paper we compute 
finite t' Hooft coupling corrections (corresponding to string theory a' corrections in the 
dual formulation) to the M = 4 plasma. Since physically equivalent results extracted 
from the equilibrium correlation functions are already available [18, 19], our analysis 
provides a test of the framework [1] beyond the supergravity approximation. We find 
that the finite coupling corrections to the shear viscosity of the M = 4 plasma extracted 
from the string theory dual to its Bjorken flow disagrees with the equilibrium correlation 
function computations. 

The paper is organized as follows. In the next section we review the results of 
[18, 19], the dual string theory model within which such computation were performed, 
and emphasize the self-consistency of the obtained results 4 . We further summarize our 
new results for shear viscosity at finite 't Hooft coupling. In section 3 we describe our 

1 For an approach alternative to [1] see [7]. 
2 Related work appeared in [14, 15, 16]. 

3 One can formulate and study questions such as dynamical hadronization in this approach. 
4 For a recent review of the Janik's framework and its clarification see [11]. 
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string theory computational framework introduced in [20]. The analysis are incredibly 
complicated and thus we focus on explaining the consistency checks on our analysis 
(fortunately, there are lots of such consistency checks). In section 4 we summarize 
equations of motion describing a' correction to the supergravity background dual to 
late-time Bjorken flow of the M = 4 plasma. We also present analytic solutions to 
these equations. In section 5 we map the results of the string theory computation 
to the near-equilibrium Muller-Israel- Stewart (MIS) theory and extract the leading 
correction to the SYM plasma shear viscosity. We conclude in section 6. Appendix 
contains some technical details; further computation details are available from the 
author upon request. 



2 Shear viscosity of the M = 4 SYM plasma at finite 't Hooft 
coupling 

Consider J\f = 4 SU(n c ) SYM plasma in the ('t Hooft) large-n c limit: we take the 
Yang-Mills coupling g\ M — > and n c — > oo in such a way that the 't Hooft coupling A, 

x = 9ym u c, (2-1) 

remains finite. We are interested in the transport properties of the plasma at strong 
coupling, i.e., A >> 1. In this limit the useful description of the Af = 4 plasma is 
in terms of holographically dual background of near-extremal black 3-brane geometry 
in the supergravity approximation of the type lib string theory. The supergravity 
approximation is exact at A = oo, in which case one finds [13, 21, 22] the speed of 
sound c s , the shear viscosity r], and the bulk viscosity £ correspondingly 

c s = ^, V = \n 2 c T\ ^ = 0. (2.2) 

In the supergravity approximation the entropy density s is 

s = ^nlT* , (2.3) 
leading to the viscosity-to-entropy ratio 

- = J- ■ (2.4) 

S 4:71 

In a hydrodynamic approximation to near-equilibrium dynamics of hot gauge theory 
plasma there are several distinct ways to extract the transport coefficients (2.2). First 
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[13], the shear viscosity can be computed from the two-point correlation function of 
the stress-energy tensor at zero spatial momentum via the Kubo formula 

77 = lim dtdx e^([T xy (x), T xy (0)}) . (2.5) 

Second [21], diffusive channel two-point retarded correlation function of the stress en- 
ergy tensor, for example, 

G tx>tx {u, q) = ~ij dtdxe^ t+ ^ z 9(t)([T tx (x), T tx (0)}) oc - , (2.6) 

have a pole at 

to = -Wq 2 , (2.7) 
where the shear diffusion constant V is 

V = ^-. (2.8) 
sT v ; 

Finally, all the transport coefficients (2.2) can be read off from the dispersion relation 
of a pole in the sound wave channel two-point retarded correlation function of the stress 
energy tensor, for example, 

G tt , tt (u, q) = -ij dtdxe-**«*'0(t)([Tu(x), T„(0)]> , (2.9) 

as 

w( , ) = e>s _ j | r 2( 1 + g. (2 ,o) 

In [18], using the Kubo formula, the finite 't Hooft coupling correction to the ratio (2.4) 
was computed. Such corrections were further extracted from the two-point retarded 
correlation functions (2.6) and (2.9) in [19]. Much like as for the infinite 't Hooft 
coupling, all three approaches led to the same result: 

7=i( 1 + f «3)A-^ + ...). (2,1) 

The computation leading to (2.11) were performed in the supergravity approxima- 
tion to type lib string theory including the leading a' correction 
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where 



d w Xv ^g 



R ~ \{94>f ~ j^y (^ 5 ) 2 + ... + 7 e-i+W + 



7 = ^(3)(a') 3 , 



(2.12) 



corresponding to 



7 = k(3) A- 3/2 



(2.13) 



in terms of the SYM 't Hooft coupling A, and 



VV — O O pmn gO h O rsfc + -o u 3 



rsp^q 



rsk ' 



(2.14) 



A somewhat uncontrollable approximation 5 used in [18, 19] (and earlier in [23] ) was 
an assumption that in a chosen scheme (2.12), the self-dual F§ form does not receive 
order (a') 3 corrections. Though the agreement of the shear viscosity extracted from 
different boundary correlation functions at order 0(7) computationally appears to be 
rather dramatic 6 , a recent work [24] suggests that such an agreement is automatically 
encoded in any gravitational model, perhaps similar to the way the thermodynamic 
relation between the energy (£), free energy (JF) and the entropy densities J- — £ — Ts 
is encoded in the horizon geometries [25]. 

In [2] the ratio of the shear viscosity to the entropy density at infinite 't Hoof 
coupling was computed for the Bjorken flow of the M = 4 SYM plasma, and the result 
was found to agree with (2.4). In this paper we compute leading order finite 't Hooft 
coupling correction to rj/s for the M = 4 plasma in the Bjorken flow. We use string 
theory effective action (2.12), which is expected to reproduce (2.11) obtained from 
the equilibrium correlation function computations. Unfortunately, we report here a 
disagreement with (2.11): 



In principle, analysis are very straightforward: 

■ we compute a' corrections to the gravitational background dual to boost invariant 
Bjorken expansion on the boundary; 

■ as in [2], we fix some of the parameters 7 of the a' corrected geometry by requiring 



time expansion; 

5 This is related with difficulties of deriving a' corrected string theory effective action in the presence 
of nontrivial RR fluxes. 

6 For example, using (2.9), (2.10) one has to take into account a' corrections to the black three 
brane temperature [23]. 

7 Amusingly, not all parameters can be fixed from the nonsingularity condition [1]. Despite this, 
our answer for the viscosity ratio (2.15) is unambiguous. 



7/ 



Bjorken flow 




(2.15) 



nonsingularity of the string theory background to order 0(t 4 / 3 ) in the late proper 
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■ we use a' corrected holographic renormalization (see [20]) to evaluate the energy 
density of the expanding M = 4 plasma; 

■ as in [11] we interpret the string theory computation in the framework of the Muller- 
Israel-Stewart theory of dissipative relativistic hydrodynamics 8 . 

Given the computational complexity of the analysis and the disagreement we claim 
to exist between equilibrium correlation function computation and the Bjorken flow 
computation of the plasma shear viscosity at finite 't Hooft coupling, the remainder of 
the paper focuses on explaining the consistency checks we performed. 

3 Computational framework 

The most difficult step in the analysis is to obtain a'-corrected supergravity equations 
governing the SYM Bjorken flow at finite 't Hooft coupling. As in [20], we first derive 
effective action for the warp factors of the gravitation background, and then deter- 
mine a'-corrected equations of motion. As explained in [11], choosing a Bjorken frame, 
fixes reparametrization of a radial and a proper time coordinates, which leads to two 
first order differential constraints (at each order in the late proper time expansion) on 
the warp factors. Not only these constraints are important consistency checks on the 
obtained system of equations, but they are also crucial to pick up correct solutions 9 . 
If one naively obtains effective action for the warp factors, the constraints are lost. 
We explain how to properly derive effective action, so that the reparametrization con- 
straints are kept. As in the supergravity approximation, the a' corrected constraints 
will provide vital consistency checks. 

3.1 Supergravity approximation: effective action, constraints and the so- 
lution 

The supergravity background holographically dual to a Bjorken flow of the Af = 4 
plasma takes form [1, 11] 

ds 2 w = g M Nd£ M di N = 

= e- 2ai - T > z) gilv (x)dx»dx v + e 6/5a(T ' 2) (dS 5 ) 2 , 

8 In [17] a correction to Bjorken expansion within MIS theory was pointed out. Such a correction, 
while modify the interpretation of the plasma relaxation time, does not affect the determination of its 
shear viscosity compare to the analysis of [11]. 

9 These constraints fix some of the integration constants. 
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dz 2 (3.2) 



where (dS 5 ) 2 is the line element for a 5-dimensional sphere with unit radius, and 
ds 2 = g av dx >x dx u = 

= 1 \ e 2a{T, Z ) dT 2 + e 2b(r,z) T 2 d 2 + e ^z) ^2 j + 

where <ix5_ = cfei + da^. The 5-form F§ takes form 10 

^5 = ^5 + ^5 , ^5 = -4 u; S 5 , (3.3) 

where u^s is the 5-sphere volume form. Moreover, the dilaton is <fi = 0(r, z). 

Notice that choosing the five- dimensional metric as in (3.2) fixes the (r, z) repa- 
rametrization invariance so that 

9zz = \ , 9rz = 0. (3.4) 

Substituting the ansatz (3.1)-(3.3) into type lib supergravity action 11 ( the 7 = 
approximation of (2.12) ) one obtains naive effective action for the scalars {a, b, c, a, 0}: 



^ef f, naive 



a(r, z), b(r, z), c(r, z), a(r, z), <j>(r, z) 



= I 



(3.5) 

7=0 



Clearly, from (3.5) one obtains only five second order PDE's, instead of seven, as in 
[11]. The reason that is happening is because the reparametrization constraints (3.4) 
are imposed directly on the action, instead on the equations of motion. 

In order the keep constraints within the effective action approach we deform the 
five dimensional metric (3.2) as follows 

ds 2 = g^dx^dx" = 

1 d? 2 

= J_ \_ e 2a(r,z) dr 2 + ^(r,*)^ 2 + ^(^21 + ^ + (1 + g fa z )) ™ . 

Z Z 1 

(3.6) 

Constraints (3.4) now correspond to 

h(T,z) = 0, g(r,z) = 0. (3.7) 



10 We normalize the five-form flux so that the asymptotic AdS radius is one. 

11 As usual, care should be taken with the self-dual five-form. The correct contribution to the 
effective action from the five-form is — j^y-Ff = — 8e~ 6a ( r ' z \ see [11]. 
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Next, we evaluate the supergravity action on (3.6), (3.3) 



'e// 



a, 6, c, a, ; ft,, 5 



(3.8) 



7=0 



While it is important to maintain full nonlinear contributions for the "physical " scalars 
{a, b, c, a, 0} in S e ff, it is sufficient to evaluate (3.8) to linear order in {h, g}. The equa- 
tions of motion derived from S e ff, followed by fixing the reparametrization invariance 
as in (3.7) 









6S, 



eff 



5a(r, z) 

SS eff 



, 



5S, 



eff 



h(T,z)=g(r,z)=0 



8a(r, z) 

6S eff 



Sb(r, z) 

SS eff 



, 



5S, 



eff 



h(r,z)=g(r,z)=0 ^ C ( T ' Z ) 



h(T,z)=g(r,z)=0 



h{T , z)=g{T , z)=0 ' ° 6<j>(r,z) 



Sh(r, z) 







h(r,z)=g(T,z)=0 
SS eff 



h( T > z )=9(T,z)=0 



5g(r,z) 



h(T,z)=g(r,z)=0 



(3.9) 



are precisely equivalent to the full set of equations derived in [11] (see eq. (3.12)-(3.18)). 

Equations (3.9) are solved as a series expansion in the late proper time r — > 00, 
but exactly in the scaling variable [1] 

(3.10) 



7-1/3 ' 

We find it useful (especially when we generalize our effective action approach to the 
full action (2.12) ) to rewrite S e ff as a functional of 

{a(r,v), b(r,v), c(t,v), a(r,v), (/>(t,v), h(r,v), g(r,v)} , 

i.e., fields depending on the scaling variable (3.10). When we do this, it is important 
to note the modification of the effective action measure 



dtdz 



t 1/3 dtdv 



Solving (3.9) in the scaling variable v and with the boundary conditions 



a(r, v), 6(r, v), c(r, v), a(r, v), 0(r, v) 



= 0, 



(3.11) 
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we find [1, 2] 

1 1 



with 



a(r, v) = Oq(v) + -^ai(v) + ^a 2 (w) + 0(r ) , 

fc(r, v) = b (v) + ^6i(«) + ^h(v) + 0(t~ 2 ) , 

1 1 

c(t,v) = c (v) + -jfiC^v) + ^c 2 (v) + C(r" 2 ) , 

a(r,v) = 0(r- 2 ), 0(r, „) = 0(t~ 2 ) , 



Q °-^ ln TT$f ' &o = c = ^ln(l + ^/3) , 
(9 — f 4 ) f 4 f 4 57 3 — f 4 



a i = ^ — ^ 5— > ci = -77 — ^ - - In 



9 -v 8 ' ' 1 3 + v 4 2 3 + v 4 ' 
3+V 

(9 + 5v 4 )v 2 ,,(9 + v 4 )v 4 „ 2 (-1053 - 171v 4 + 9v 8 + 7v 12 )v 4 
°2 —rw77\ — -*c - 6 ^77, — "ST + V 



b i = ~3v , , 4 - 2ci , 



12(9 -w 8 ) 72(9 -v 8 ) ' 6(9 -v 

1 , v^-v 2 3^ 2l 3-v 4 
H — In — — 77 In , 

8V3 V3 + v 2 4' 3 + v 4 ' 



8^2 



7T 



2 v 2 (9 + v 4 ) ^ v 4 „ 2 (-9 + 54v 4 + 7v> 4 



288^ 12(9 -v 8 ) 72(3 + v 4 ) 1 6(3 + v 4 )(9 - v 8 ) 
1 , V3- v 2 1 ,„ 3-v 4 



H ^ + — (C + 66?f)ln o 

8^ v^ + v 2 72 v /; 3 + v 4 

. , ... ,....,...,» / r,/^-,- 2 

+ — — I In — r In tttt- — — li 2 



24VS\ VS + v 2 4(3 + v 4 ) 2 z y (v^ + v 2 ) 2 

(39 + 7v 4 )v 4 1 , V3- 



2 4 

b 2 = - 2c 2 + — - T\ + C— - + - (f - + — 7= In 



4(3 + v 4 ) 24(3 + v 4 ) ' 2(3 + v 4 ) 2 8^ 



3. 2l 3-v 4 
+ l 71 ln 3+V' 
for arbitrary parameters {fj, C}. 
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3.2 Beyond the supergravity approximation 

To go beyond the supergravity approximation we deform proper late time expansion 
for the physical fields (3.12) as 



a 2 (v) + 762 (t,v) , 



a(r, 


«) = 


b(r, 


«) = 


c(r, 


v) = 


a(r 


v) = 


0(r 


v) = 



r 2/3 
1 

^273 
1 



-2/3 



r 4/3 
1 

TITS 
1 

7473 



1 



1 



T 



4/3 



r 2/3 f^V'>" 

^273 7Mt,v) + ^702 (r,u), 



where we assumed that as A — > 00 

aj(Ar, i>) -> aj(r, i>) , 6j(Ar, u) -> &j(r, t>) , q(At, d) -> Cj(r, u) , 
di(Ar, u) -> dj(r, u) , 0j(Ar, u) -> 0;(r, u) , i = 0, 1, 2 . 

Substituting (3.16) into (2.12) we can obtain effective action for 

{ai,bi,Ci,ai,(f)i} . 



(3.16) 



(3.17) 



(3.18) 



We need to expand the supergravity part of (2.12) to quadratic order in the fields 
(3.18), while it is sufficient to evaluate the 'je~ 3 ^ 2 ^W term in (2.12) to linear order 
in (3.18). As explained in the previous section, in order to obtain reparametrization 
constraints we further need to evaluate the whole action (2.12) to linear order in {h, g}. 
Resulting effective action can be organized as follows 



C7 



eff 



di,bi,Ci,ai,^i; h,g;r),C 

I ^[2] + 7^[l] 



dr I dv Clff 



[i] 



di, b h di, a h fa; h,g;rj,C 

^2 £ -5/3 £-5/3 
I *~[2] + '^[l] 



(3.19) 



where, given (3.17), the superscript in C indicates the leading scaling behaviour as 
r — > 00, i.e., 

C q = 0(r q ). (3.20) 
Notice that we suppressed explicit dependence on {rj, C} in the decomposition of Clff- 
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In the rest of this subsection we present further decomposition of C q , and explain 
the structure of the consistency checks on the example of the equation of motion for 

In the following expressions all £'s do not scale as r — > oo — the explicit scaling 
dependence is factored out. A superscript indicates whether a term comes from the 
supergravity part of (2.12), or from the (9(7)-correction in (2.12); the subscripts indi- 
cate order of fields on which a lagrangian term depends. For example, the subscript in 
£(2,1) implies that £(2,1) is strictly bilinear in order order-2 and order-1 fields: 



'(2,1) = £(2,1) 



a 2 , h, c 2 , «2, 02 ) , (a>i, h, ci, di, <p^j 



while the subscript in £(1,1) implies that such a term is strictly quadratic in order-1 
fields 

ai,6i,ci,&i,0i) , (ai,&i,ci,di,0 



£(1,1) = £(1,1) 

A single subscript indicates the order of fields on which the corresponding lagrangian 
term depends linearly. Finally, no subscript indicates that the lagrangian terms are 
independent of (3.18), but are linear in either h or g. Notice that lagrangian terms with 
a single subscript coming from the supergravity part of (2.12) are necessarily linear 12 
in either h or g. The latter implies that reparametrization constraints at leading order 
in 7 appear as 0(7), and thus in the lagrangian terms which are either bilinear or 
quadratic in the fields (3.18) we can set h — g — 0. For the same reason we can set 
h = g = in the lagrangian terms of order 0{^ 2 ) coming from the C(7)-correction in 
(2.12). 

3.2.1 Order- 
We find 



,-1/3 _ 1 



12] - 1/3 



-1/3] r SUGRA , [-1/3] r W \ 
^(0,0) + ^(0) J » 



, f (3-21) 

^-1/3 _ 1 / [-1/3] rSUGRA , [-1/3] r W \ 
^-[1] ~ r l/3 I ^(0) "T J ' 



12 This is simply a check (which we verified to be true) that the supergravity background (3.13)-(3.14) 
solves supergravity equations of motion. 
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3.2.2 Order- 1 
We find 



(3.22) 



r -3/3 _ 1 / [-3/3] rSUGRA , [-3/3] A 
^[2] ~ r ^ ^(0,1) + ^(1) y 

i ! ( [-5/3] r SUGRA , [-5/3] r SUGRA , [-5/3] A 
+ 7B73 I MO,!) + ^(1,1) + ^(1) I , 

£-3/3 = 1 ^ [-3/3] £ S!7GiM + [-3/3] £ Sl7GflA + [-3/3] (3 33) 



5. £5 Order- 2 
We find 



/*-5/3 / [—5/3] rSUGRA . [-5/3]^ ] 

^[2] - r 5/3 I ^(0>2) + ^(2) J 

1 ! / [-7/3] rSUGRA , [-7/3] rSUGRA , [-7/3] /*VF \ 

+ r 7/3 I ^(0,2) + ^(1,2) + ^(2) I 

1 ! I [-9/3] rSUGRA , [-9/3] rSUGRA , [-9/3] rSUGRA , [-9/3] /-IV A 

"^9/3 I ^(0,2) + ^(1,2) + ^(2,2) + ^(2) J 



(2) ; > 

(3.24) 



/>-5/3 _ 1 / [—5/3] rSUGRA . [—5/3] rSUGRA , [-5/3] rSUGRA , [-5/3] /*WM /o ocr\ 
^[1] - 7B73 ^ ^(0) + ^(1) + ^(2) + ^ J ■ 

3.2.^ Equations of motion for 0,2 and consistency check on lower order solution 
Consider a variation 

Given (3.19), (3.24), we have 13 

4 ( [ - 7/3] X(o) + [ - 7/3] X(i) + h7/3] x) (3.27) 

-^3 ( [ " 9/31 J (o) + h9A %) + h9/3l ^( 2 ) + h9/31 ^) > 
13 The term [-?>/3\£SUGRA j oes not con t r ibute as it can only be bilinear in both S2 and either h or 



+ r 7/3 



9- 
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with obvious notations. For example: 

^ M/3,I <»> = ^b)J dt J dv l " 7/3|£ « M } ■ (3 - 28 > 

Thus, a single variation (3.26) not only will produce the equation of motion for a 2 , 
namely 

0= ^% 0) + [ ^ 9/3l J (1) + [ - 9/3] T( 2) + [ - 9/3] J, (3.29) 
but also would produce a consistency tests on the order-0 and order-1 solutions: 

=HV3]J + [-5/3] J 

( j (3.30) 
= I- 7 / 3 ! J (0) + [- 7 / 3 ]j {1) + I- 7 ^ J, 

correspondingly. It is easy to understand the origin of consistency constraints (3.30): 
a field deformation at order-2 can be considered as a variation of a field deformation 
at order-0 or order-1; thus, once equations of motion at lower orders are solved, such 
variations must vanish. 

Altogether, there are 10 consistency constraints when deriving equations for order-2 
field deformations, and 5 consistency constraints when deriving equations for order-1 
field deformations. We explicitly verified that all such constraints are satisfied. 



4 Equations of motion and solutions for (3.18) 

We obtain equations of motion (including the constraints) at each order % = 0, 1, 2 
for (3.18) as explained in the previous section. It is important that in deriving these 
equations of motion a field is assumed to depend both on {r, v}. Once the equations 
of motion are computed, we can assume 

{a h k, ti, a u 0j}(r, v) = {a h k, q, a h 4>i}{v) . 

All the equations must be solved with the boundary conditions 



a,i(v), bi(v), Ci(v), &i(v), 4>i(v) 



= 0. (4.1) 

v-<-0 



4.1 Order-0 



To leading order we find the following system of equations 

1 tn 5v 4 - 9 v 5v 4 - 9 f . 933120(4w 8 - 33w 4 + 36)v 10 

13 



.„ l.„ 5v 4 + 9 5v 8 + 27„ 311040v 10 (> 8 -45t; 4 + 36) , 4 nS 

°= c » + 2°« + 2F^ a ° + ^i) Co ■ (4 - 3) 



.„ t„ 5w 4 + 9 5v 8 + 27 f/ 5v 8 + 27„, 

=cg + a'J + bl + — — -a' + g b' + 8 c 
(f 4 — 3)v v(v s — 9) v (v 8 — 9) 

622080w 10 (4t; 8 - 45w 4 + 36) 

(3 + v 4 ) 8 ' 

U 3 (t; 8 -Hh; 4 -3) A/ _ 3(5i; 4 -3U _ 3(^ 4 - 3) 
U_C ° + 2 &0+ ^ 8 -9) C ° v ( v 8-9) 2(3 + V > a ° 
3732480t; lt V - 3) 2 
+ (3+^j 8 ' 

1- (v 4 - 3) 2 ^ 

-c + -6b + 2{v4 + 2v2 + 3){v4 _ 2v2 + 3 fo 



(4.4) 



(4.5) 



2799360(w 4 - 3)v 
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(4.6) 



(v 4 - 2v 2 + 3){v 4 + 2v 2 + 3) (3 + w 4 ) 7 ' 



5v 8 + 27„, 32 „ 2332800w 14 
f (v s — 9) ir (3 + f 4 ) 8 

7 „ 5w 8 + 27 ? . 5598720w 14 

°=« + ^-W (4 ' 8) 

Although the system (4.2)-(4.8) is overdetermined (we have 7 ODE's for the 5 func- 
tions), it is straightforward to verify that it is consistent. 

Most general solution to (4.2)-(4.8), subject to the boundary conditions (4.1), is 
parametrized by three arbitrary integration constants 14 {8i, f3i, fa}: 

72576 585144 216 1714608 864 1714608 
a ° ~ (3+^F ~ (3 + v 4 ) 4 ~ 3+^ + (3 + t; 4 ) 5 ~ (3 + v 4 ) 2 ~ (3 + v 4 f g) 

t 223560 25920 664848 664848 5 Y 1 . 

+ 7^^-7^7—^-77^7^ + 24+^!, (4.10) 



u (3 + v 4 ) 4 (3 + t; 4 ) 3 (3 + v 4 ) 6 (3 + t; 4 ) 5 4(3 + v 4 ) 12 

c = k, (4.11) 
. ( ( v 4 1 \ 3 - v 4 1 \ 25 / 4 9 \ , 3 - v 4 75 

^ + ^J + ^r + ^) ln 3T^ + Y 

225 675 22275 36450 36450 ^ ' ^ 

+ 2(3 + v 4 ) 2 + (3 + w 4 ) 3 ~~ 2(3 + v 4 ) 4 + (3 + v 4 f ~ (3 + v 4 ) 6 ' 



14 Altogether we expect 10 integration constants; setting nonnormalizable components of the fields 
to zero fixes 5 integration constants; the two constraints fix another 2 integration constants. 
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f(3 2 45 \ , 135 405 810 25515 



v 24 4 / 3 + w 4 2(3 + w 4 ) 2(3 + w 4 ) 2 (3 + v 4 ) 3 (3 + w 4 ) 4 , . 

87480 87480 1 ' ' 

~ (3 + w 4 ) 5 + (3 + w 4 ) 6 ' 
For a generic choice of /3i the curvature invariants of the ten-dimensional geometry 
(3.1) will be singular as v — > 3^ 4 . For example, 

K^pxRr* = ~1 (1350 + A) ln(3 - t; 4 ) + finite . (4.14) 
However, once we chose 

p 1 = -1350, (4.15) 

the curvature invariant (4.14) is finite. We verified that for (4.15) all other curvature 
invariants (including those discussed in [11]) are finite. 

Notice that unlike studies of the string theory dual of the Bjorken flow in the su- 
pergravity approximation [1, 2, 3, 11], here the background geometry is not completely 
fixed from the nonsingularity condition — two arbitrary integration constants {Si, fa} 
remain. Though (3 2 does not affect the holographic stress energy tensor of the back- 
ground geometry, the latter does depend on Si. Additional arbitrariness will appear at 
higher orders, however, the string theory computation when properly matched to the 
near-equilibrium SYM dynamics provides unambiguous predictions for plasma trans- 
port properties. We further comment on the physical origin of the arbitrariness in 
section 5. 

4.2 Order-1 

Using results of the previous subsection, including (4.15), we find the following set of 
equations for the next-to-leading order in the late proper time expansion at order 0(j) 

^ + ¥ + + |^!f + 5 <"> • < 4 - 17 > 

= c ' + a ' + b'[ + 4 a + 8 6 + 8 c + 5 (1>3) , (4. 18) 

(f 4 — 3)v v(v s — 9) v (v s — 9) 



,2 



I?,, (v 4 -15> 3 ~, 5w 8 -30t; 4 -27^ 3(w 4 - 3) ^ 24w 2 
° " Cl + 2 1 + V 8_ 9 6 i + v(u 8 _ 9) c i " 2 (3 + t;> ai " ^^9 1 
48t; 2 A 

-Ci + 0(1,4) , 



(4.19) 



v 8 -9 



15 



0= ^ + 2(^- J+^ + 2^ + 3) a ' + ^ + ^" (420) 
5tj 8 + 27., 32 „ , 4MS 

O=0'/ + |^0' 1+ ^ 7)! (4.22) 

where the source terms {«S(i,i), • • • ,5(1,7)} are given in Appendix A. As before, while 
the system (4.16)-(4.22) is overdetermined, we explicitly verified that it is consistent. 

Most general solution to (4.16)-(4.22), subject to the boundary conditions (4.1), is 
parametrized by three new arbitrary constants {5 2 , fa, Pa}' 

72013536?} 61725888?} 31189536?} 3r)(288 + <&i) 59 8 75207) 
ai ~~ (3 + ^ 4 ) 6 + (3 + *T + (3 + ^)« + (v*-3Y ' JsTW , . 
3?}(1728 + gi) 576?} + f £ 2 |^ - 576?} 445824?) 1 ^ 
~ 2(3 + ?; 4 ) 2 + 2(?; 4 - 3) + 2(3 + v±) + (3 + ~ Q V 1 + 18 ' 

? /^ , Mi 3-u 4 86477 - ^ + 1 144t) + ±7^ 3024?} 

&1 = - 12077 - -77^1 + — In — — 2_ + 



V""'' 6 ' 18 J " 3 + t; 4 3 + t; 4 v 4 - 3 (3 + tj 4 ) 2 

37}5i 1168 73287} 1071367} _ 23934528?} 2091744?} 27597024?} 
+ 2(3 + t; 4 ) 2 ~ (3 + t; 4 ) 5 ~(3 + t; 4 ) 3 ~ (3 + t; 4 ) 7 + (3 + tj 4 ) 4 + (3 + t; 4 ) 6 

(4.24) 



1 ~A a. 52 
-6^ 1+ 18' 



2i = (60?} + | - 1^) In. 



. . <r 72?} + jgft 432?} - | - jgft 1 

179712?} 23934528?} 2336688?} 28086912?} (1512 + §<5i) fj 
~(3 + v*) 3 ~ (3 + v 4 ) 7 + (3 + t; 4 ) 4 + (3 + w 4 ) 6 + (3 + t; 4 ) 2 

12177216?} 5 2 
~ (3 + ?; 4 ) 5 + 18 ' 

n ( ( v 4 1 \ , 3 - v 4 1\ 25?} / 4 9\ 3-« 4 75?} 
ai = -^U^ + 96^J ln 3T^ + lliJ + — + ^J ln 3T^ + — 
450?}?; 8 (?; 4 -3)(?j 8 + 24?; 4 + 9) 
+ (3+W ' 



(4.26) 



/3 4 /3 2 r}\ , 3 - v 4 r)(3 2 W 2 1487160?} 3149280?} 

In ■ 



24 12 y 3 + ?j 4 4(3 + v 4 ) 4(w 4 -3) (3 + v 4 f (3 + w 4 ) 7 
405?} 135?} 135?} 2430?} 218700?} , 3674160?} 



(3 + ?j 4 ) 2 2(3 + w 4 ) 2(w 4 -3) (3 + ?j 4 ) 3 (3 + w 4 ) 4 (3 + v 4 f ' 

(4.27) 
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For a generic choice of /3 3 the curvature invariants of the ten-dimensional geometry 
(3.1) will be singular as v — > 3 1 / 4 . For example, 

n^xTZ^ = y 7 (& - 5400r}) In (3 - t; 4 ) + finite . (4.28) 
However, once we chose 

p 3 = 5400r} . (4.29) 

the curvature invariant (4.28) is finite. We verified that for (4.29) all other curvature 
invariants (including those discussed in [11]) are finite. 

At this stage we have four arbitrary integration constants: {5 1 , /3 2 , #3, At}- 

4.3 Order-2 

Using results of the previous two subsections, including (4.15) and (4.29), we find the 
following set of equations for the next-to-next-to-leading order in the late proper time 
expansion at order 0{j) 

=« + + WT^fi + wrwfi + Sftl) • (430) 

^ + 2 ai + + S^3)° 2 + 5<2 ' 2 » ' <4 ' 31) 



9 + 5w 4 5v 8 + 27~, 5v 8 + 27 

t;K-3) a2 + ^ _ 9 ) 2 + ^8_ 9 ) r 



=% + a'i + + J A ^ a' 2 + — - b' 2 + ^ c' 2 + 5 (2 ,3) , (4.32) 



n ,„ U„ 3 - 3 2v 8 - 15v 4 - 9 f/ 7w 8 - 30v 4 - 45 48v 2 ? 
=Co + -b~ y -^-aL H — r — b' H — c' 9 b 2 

2 2 2 2(3 + v 4 )v 2 v(v 8 -9) 2 v(w 8 -9) 2 v 8 -9 

96v 2 A 
~ ^9 C2 + lS ( 2 ' 4 ) ' 

(4.33) 

1- (V - 3) 2 

=c' 2 + -6o + — - 1 7T \~r~i 5 ra 2 + S( 2 5) , (4.34) 

2 2 2 2(w 4 -2w 2 + 3)(t; 4 + 2t; 2 + 3) 2 ( ,) ' V ; 

5^8 _|_ 27 32 

=d 2 + Tl - — «2 + 5(2,6) , (4.35) 

f (f 8 — 9) f 2 

EU;8 _|_ 27 - 

°^2 + ^9)0 2 + «S (2 , 7 ), (4.36) 

where the source terms {5( 2 ,i) •• • ,5(2,7)} are given in Appendix B. As in previous 
subsections, while the system (4.30)-(4.36) is overdetermined, we explicitly verified 
that it is consistent. 
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Solving (4.30)-(4.36) is quite complicated. Fortunately, we do not need a complete 
solution. Our ultimate goal is to determine fj from the nonsingularity of the ten di- 
mensional metric curvature invariants to order 0(7) and to order (9(r~ 4//3 ) in the late 
proper time expansion. Thus we evaluate metric invariants first, find what field com- 
binations affect the singularity as v — > 3 1 / 4 , and then solve just for those combinations 
of field. 

Up to now, all our equations and solutions are exact in fj. We assume now that 

fj = fjo + 7771 + 0( 7 2 ) , (4.37) 

and evaluate background curvature invariants to order 0(7) near 

x = 3 1/4 -v. (4.38) 

We use explicit solutions (3.15), (4.9)-(4.13) and (4.23)-(4.27) to evaluate curvature 
invariants. At this stage we use equations of motion (4.30)-(4.36) to eliminate the 
derivatives (if possible) of {a 2 , b 2 , c 2 , a 2 } from the curvature invariants. 



4.3.1 K 

For Ricci scalar we find (there is no dependence on order-2 fields here) 

3 l&j* -y/3 (\ 2 1 13^ 1 3V 4 5 \\ n 
U = -2 1 —^ {-*-W-*-—-* + —x-) + fimt6 ' X -° + - 

(4.39) 

From (4.39) we find that Ricci scalar of the string theory geometry is nonsingular as 
x — > 0+ when 

^ = ^571 > ( 44 °) 
which is precisely the condition found from the nonsingularity of Riemann tensor 
squared [2], as well as higher curvature invariants [11], in the supergravity approxi- 
mation to the string theory dual of the M = 4 SYM Bjorken flow. The difference here 
(compare to [2, 11]) is that 770 is already fixed by requiring the nonsingularity of the 
Ricci scalar 15 . 

15 Notice that Ricci scalar vanishes in the supergravity approximation [11]. 
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4.3.2 n, upX n^ x 

A bit more work is necessary to determine the nonsingularity condition of Riemann 
tensor squared at order 0(7). Generalizing the notation of [11] 

= {lt ]SUGR \v) + 7 2?%)) + ^ (l? SUGRA (v) + 7 Xf '%)) (441) 

+ ^ (4 2]5C/Gi %) + 7 4 2] ») + 0(-- 2 ) + 0(f) . 

With (4.40) and the lower order solutions, all terms in the decomposition of in 
(4.41) except for I^ W (v) are finite as v — > 3^ 4 . We find 

. [2 ]vf _ 2304w 5 (t; 4 - 3) f , ( ^ 192(5v 16 + 60w 12 + 54w 8 + 540w 4 + 405) 



= ' (3 + ^)3 ^ (3?W ^ 



+ 5X 2 ( V ), 



(4.42) 



where 

h{v) = c 2 (v) + ±b 2 (v) , (4.43) 
and 5T 2 (v) does not depend on order-2 fields, and can be evaluated explicitly 



Vs 1 

6 x 2 

(4.44) 



Sl2 = (_ 399 V3 + ^ - + Wl ) x ( 1 - JL L 3 _ 

3 1 / 4 (72^ 3 3 / 4 771 - 4v/2 3 3 / 4 5 2 + 39 ^ - 270) 1 \ 

H 7^ f= ; f= ; — s + finite , 

2 (-7182 + 15 5 ± - 4V2 3 3 / 4 5 2 + 72 ^ 3 3 / 4 m) x J 

as x — > 0+. 

From (4.30) one can obtain a decoupled equation for / 2 as defined in (4.43). The 
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resulting equation is possible to explicitly solve for f 2 : 

j . . 2412600422477 V 13961808?/, 9 9 ^, 34992?/, „ 

*M = (* + 3)» + 1^T3F (2 " + 1932,} - c) - FTf 1 * 

243?/ 

+ 320520?f - 399C) + - (28221?/ + 8578560?f - 20600C) 

81 ?; 3 97?; 3 

; (20823?/ + 3918720?? 2 - 185 60C) + — rr(7561u 2 + 485 1 7? 2 



2(?/ + 3) 5V ' ' ' 8(?/ + 3) 

3 

+ 921600rf - 11520C) — ^ (-26307?/ + 8u 2 <?i - A5 X C + 407808?} 2 

16(f 4 + 3/ 

, P , 2xx / 35831808?/ - 1990656?/ 

-1327104?/ + 29859840?) 2 (8709120?) 2 - 290304?/) 

+ (?/ - 3) 4 + (?/ - 3) 3 

If 3 
+ 7^ — 77? -23472?/ + 1029024?f - -v% - 275 1 fj 2 
(?j 4 — 3) z \ 2 

+ (Vf - It; 2 *! + 25488rf - 6f)5 2 + ^?/ 

10425^ ?/ r \ 

+ ^e^ arctan vl + 7' 

(4.45) 

where 53 is a new arbitrary integration constant. The solution (4.45) is exact in rj. 

Unfortunately, we can not present an explicit solution for d 2 (v). We can prove 
though that a 2 (v) can be chosen to be finite as v — > 3 1 / 4 , while having a vanishing 
nonnormalizable mode as v — > + . Indeed, the most general inhomogeneous solution 
of (4.36) as rr — > 0+ takes form 

a 2 (x) = Co + d In +* + ^ 3 3 / 4 C + g| 3 1 / 4 ) + O(o; 2 lux) , (4.46) 

where Co and Ci are the two arbitrary integration constants. On the other hand, the 
homogeneous solution to (4.36) is 

^hom( n ,\ _ rhom V 8 + 9 | nhom ( ( V 4 1 \ , 3 — V 4 1 



") = C m — r- +Cr m 7^ + 7^ ^7^—7 + 777 - (4-47) 



?/ 1 V.V.864 96?// 3 + ?/ 144 

leading to asymptotic expansion as v — > 0+ 



i nhom 

a h 2 om ( v ) = 9 — + ^ 4 - 7^777 + 0(O , (4-48) 
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and 

a h 2 om (x) = 6C feom +^Cf om ln2-^Cf om ln3+^Cf om +^C{ lom lna;+C'(x) , (4.49) 
as x — > 0+. Thus, given (4.46), we can take 

C\ om = -144Ci, (4.50) 

and tune 16 C^ om so that 

a v 2 + a h 2 ° m (v) 

both have only a normalizable mode as v — * 0+, and be finite (along with its derivatives) 
as x — > + 

We now have all the necessary ingredients to determine 171 from the nonsingularity 
of l [ 2 ]W (v). We find that if" 

{y) is finite as v — > 3_ , provided 

7)1 = ^3lT (7182 " 15 5l + ^ ^ • (451) 
Notice that while 770 (4.40) is determined unambiguously from the nonsingularity condi- 
tion of the background geometry, the absence of singularities is not a powerful enough 
constraint to fix fji. As we already mentioned, this fact will not preclude us from 
computing a definite ratio of shear viscosity to the entropy density. 

4.3.3 n^w* 

Analysis of the square of the Ricci tensor can be performed in the same way as for the 
Riemann tensor square. We find 

T> W fl ■+ 1920 "2(u) Q l/4 u 

TZ^JRr = finite - 7 —q^ , v -> 3 J , (4.52) 



T 



were we explicitly indicated dependence on order-2 fields. We argued above that &2 (v) 
can be chosen to be finite as v — > 3 1 / 4 ; this would guarantee the nonsingularity of 
n^JRr to orders 0{j-^) and 0(7) . 

^.5.^ Higher order curvature invariants 
As in [11] we denote 

' ^i/pA = ' nivijxv ■ ' ^ 1 pA 5 (4.53) 
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This is always possible since £(2.6) = C( w8 )> see (B.9). 
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where 

^ [0] , vp x = n, vpX . (4.54) 
We further define higher curvature invariants l^ n \ generalizing (4.41): 



= (4 2 " ,Sl ' GM («.) + 7 lP W (v)) 4- +, (IT"!.) + 7 O)) (4.55) 
+ 7^73 (lP SUCRA (v) + 7 if '%)) + 0(r" 2 ) + 0(7 2 ) ■ 

Given complexity of analysis, we checked at order (9(r~ 4 / 3 ) only 17 nonsingularity 
of JW. Keeping explicit only fields at order-2, we find 

4 4] %) = - 552 ^ + 5( ; 4 4 )7 " 3) («» - ^ + - 36, 4 + 8l) 

663552^(^-3)3 _ 1536 / 

(3 + v 4 ) 7 21 ^ (3 + v^\ (4.56) 

+ 33480v 20 - 72738v 16 + 301320v 12 + 556956v 8 + 87480v 4 

+ 32805 ) & 2 (y) + SX A {v) , 



No new analysis of order-2 fields are necessary here to verified that (4.56) is nonsingular 
precisely for (4.51). 

5 Shear viscosity for the Bjorken flow of J\f = 4 SYM plasma 

In previous section we analytically evaluated a'-corrected supergravity background 
dual to a Bjorken flow of M = 4 SYM plasma at finite coupling to order 0(t~ 2 ^ 3 ) in 
the late proper time expansion. We can now extract the boundary energy density e(r) 
from the one-point correlation function of the boundary stress energy tensor using the 
(^'-corrected holographic renormalization developed in [20]. We confirmed that the final 
expression for the energy density can be evaluated as in the supergravity approximation 
[1, 2, 11]: 

N 2 2a(v,r) 
e(Tj = — — — - hm — — -j— . (5.1) 



17 We studied more general higher curvature invariants at lower orders in the late proper time 
expansion. 
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Using (3.16), (4.9), (4.23), (4.37), (4.40) and (4.51) we find 

iV 2 (6 + 576 7 + 7 ^) 1 iV 2 2 1 / 2 3 1 / 4 (15667 + 8 + 7^1) 1 ^ . _ 8/3 , 



12tt 2 r 4 / 3 48tt 2 



T 



(5.2) 

Notice that even though r\\ depends on an arbitrary integration constant 62, see (4.51), 
such a dependence disappears in (5.2). 

The string theory result (5.2) should now be interpreted within Muller-Israel- 
Stewart theory [8, 9]. For the Bjorken flow of the Af = 4 SYM plasma we expect 

[11] 

e9 auge^ = 3 ^2^4 ^ + ^ J_ _ ^ rf* \* A (1 + 15 7 ) \ + 0^^) , (5.3) 

8 r 4 / 3 r 2 

where A is an arbitrary scale, related to the initial energy density of the expanding 
plasma, and 

A=-, (5.4) 
s 

is the ratio of shear viscosity to entropy density. Notice that in (5.3) we have to use 
Af = 4 equation of state at finite 't Hooft coupling 18 , see [23]. Matching the leading 
terms in the late time expansion in (5.2) and (5.3) we find 

1 



A — ^1 + 8I7 + -7 *) + 0(Y) , (5.5) 
which implies that 

A = i- (l + 120 7) +0(7 2 ), (5.6) 

leading to the result quoted in (2.15). Notice that the remaining arbitrary integration 
constant Si got "absorbed" into a definition of an arbitrary scale A. The physical 
observable, i.e., the ratio of shear viscosity to entropy density (5.6) is evaluated unam- 
biguously. 



6 Conclusion 

In this paper we evaluated finite 't Hooft coupling correction to strongly coupled Af = 4 
SYM plasma undergoing 1+1-dimensional Bjorken expansion. Such finite coupling 
correction corresponds to string theory a' corrections in the supergravity dual to the 
boost invariant plasma expansion proposed by Janik and collaborators [1, 2, 3]. We 

18 We would like to thank Romuald Janik for pointing this out. 
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extracted the ratio of shear viscosity to entropy density at large, but finite, 't Hooft 
coupling. Our result is in conflict with analysis of shear viscosity of strongly coupled 
SYM plasma extracted from the equilibrium two-point correlation functions of the 
boundary stress energy tensor [18, 19]. 

Though we do not currently understand the source of the discrepancy, there are 
two potential outcomes. 

■ First, the nonsingularity approach of [1, 2, 3] is simply inconsistent — it was shown 
in [11] that it is definitely inconsistent within the supergravity approximation at high 
orders in the late proper time expansion ( though it appears to be well-defined in 
the first two orders necessary to compute shear viscosity of Af = 4 plasma in the 
supergravity approximation ). 

■ The discrepancy might be due to using incorrect effective string theory action (2.12), 
as a dual to Af = 4 dynamics at large, but finite, 't Hooft coupling. Specifically, the full 
set of a' corrections to RR fluxes is unknown. It is conceivable that proper inclusion of 
a' corrections to the RR fluxes would reconcile two results (as it must be). A starting 
point for such analysis could be the a' effective action proposed in [26]. 

We believe resolution of the puzzle presented here deserves further study. 
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A Appendix: Source terms for eqs. (4.16)-(4.22) 



72fjv 6 5 1 1119744(20v 16 -34(h; 12 + lOlOv 8 - 555v 4 - 9)r?v 6 
' S(1 < 1) -"(3TW + {3 + ^ (A - 1} 



72(v 16 + 54v 8 + 72v 4 + 81)f)v% 373248r?v 6 A 100in i 6 

<~>(1 21 — / — 5 7^~a 77! ,, n / a 777T x 113319v 

(,} (v 8 -9) 4 (3 + v 4 ) 9 (v 4 -3) 4 V 

- 55647v 4 + 159732v 8 - 180090v 12 - 37251v 20 + 2187 - 532v 28 + 12v 32 



+ 6750V 24 



(A.2) 
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144(V 6 + Qv 12 + 5Av 8 + 18v 4 + 81)r^ 6 <5i 746496r^ 6 / 32 

5(1 ' 3) " 0^Tg)4 (3 + ^)V-3) 4 X V 

- 604v 28 + 6840w 24 - 37017w 20 + 11375k; 16 - 179118v 12 + 139806v 8 



- 30861 r J ■+ 2187 . 

(A.3) 



12^ 4 + 15X^ + 9)^ 186624^ 6 / 4 

+ 15663w 16 - 50778w 12 + 74304w 8 - 37098v 4 - 1215^ 



(A.4) 



24v 7 f)S 



i 



5(1,5) (v 8 - 9) (v 4 - 2v 2 + 3) (v 4 + 2v 2 + 3) 

373248 ^ x(V -52(kv 16 (A.5) 



(v 4 + 2v 2 + 3)(v 4 - 2v 2 + 3)(v 4 - 3) (3 + w 4 ) 8 
+ 1455w 12 - 1269w 8 - 189w 4 - 81 



_ 259200(7k 8 + 639 - AhQv A )v lA i) 

5(1 ' 6) = (3 + V *)V-3) ( } 

_288^ 10 /? 2 44789760w 14 r/(t; 8 -9t; 4 + 27)(t; 8 -3t; 4 + 3) 
(1,7) ~ (v 8 - 9) 3 (v 4 -3) 3 (3 + t; 4 ) 9 ( ' ' 

B Appendix: Source terms for eqs. (4.30)-(4.36) 

We find it convenient to express the source terms at this order implicitly in the back- 
ground warp factors {a2,&2,C2}, while using equations of motion for the supergravity 
background to eliminate higher than second derivatives 19 . For convenience, we present 
these equations here [2, 11] 



„ _ 9 + 5v 4 , lQv 3 32v 3 , 4(2v 4 + 9)v 4 



v(v 4 -3) 1 v 8 -9 z v 8 - 9 2 3(v 4 - 3)(v 8 - 9) 
A32v 6 fj 2 (v w + 2Av 12 + 30v 8 + 216v 4 + 81) 
(v 8 - 9) 4 ' 

5u 8 -8u 4 + 27, 16v 3 , A32v 6 fi 2 (v 12 + llv 8 + 63v 4 -27) 
2 ~ v(y 8 - 9) 2 + ^9° 2 (v 8 -9) 3 (3 + v*) 

4(2t; 8 - 9v 4 - 63)?; 4 
3{v 4 - 3) 2 {v s - 9) ' 



(B.1) 



(B.2) 
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The highest derivative is the fourth. 
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„ _ 8v 3 5v 8 - 16v 4 + 27 , 432v 6 rf(v 12 - 13v 8 - 9v 4 - 27) 

° 2 ~v^9 2 v(v 8 - 9) ° 2 (3 + v 4 )(v 8 - 9) 3 

4(v 8 - 9v 4 - 18)v 



4 1C ^,4 (B.3) 



3(v 4 -3) 2 (v 8 -9) ' 
Of course, (3.15) solves (B.1)-(B.3). 



_4320v 7 (119v 16 - 837v 12 + 4776v 8 - 7533v 4 + 9639) , / 2v 3 5 1 
S{2 ' 1] ~ (v 4 - 3) (3 + v 4 ) 7 ° 2 + V(s + ^ 4 ) 2 

864v 3 

- — j-y-^ — -^(125v 28 - 15725v 24 + 145095v 20 - 669570v 16 

\o -\- v ) (v 3) 

+ 1303911V 12 - 1275669v 8 + 96957v 4 + 8748)^) (b' 2 + 2c' 2 ) 

/ (2v 8 + 3v 4 + 9)v 4 432v 6 (2v 12 -21v 8 + 36v 4 -81)r/ 2 \ 
+ V (^ 4 -3) 3 (3 + v 4 ) 2 + (3 + v 4 ) 5 (v 4 -3) 3 ) 1 

" 72 K~ 2 C +9y I 2 + .^w^uo >< f-87058152,- 
(v 4 - 3) 2 (3 + v 4 ) 4 (v 4 -3) 6 (3 + v 4 ) 10 V 

- 22973922v 24 + 61231383V 20 - 39438171v 8 + 92195901v 12 + 7794468^ 
+ 177147 + 6804783w 28 - 1074783v 32 + 126421v 36 + 132v 44 - 6007v 40 J 



_ ^! x f v 40 + 2386v 36 + 5073v 32 + 50304v 28 + 2354787v 24 

(3 + v 4 ) 7 (v 4 - 3) 6 V 

- 3096360v 20 + 26563545v 16 - 11763144v 12 + 26305236v 8 - 118098v 4 

- 118098^ 

(B.4) 
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2v z 5i 288v'- 



- [j^-i + (HTrFf ' 535 '" " 5190l ' 2 ° + 34755 " 16 " 929251,12 

\ 288v 7 / 

+ 160542v 8 - 109593v 4 + 8748) )a' 2 + — -^—-^ — -- 473v 24 - 18797v 20 

/ (3 + v ) (v — 3) \ 

+ 177291v 16 - 787422v 12 + 1595619v 8 - 1522557v 4 + 344817^ 6' 2 

+ (- , 2 t V \ + 7 Tf, ,^ 24 ~ 14909, 20 + 162387v 16 

V (v 8 -9) 2 (3 + v 4 ) 7 (v 4 -3) sV 

- 751134v 12 + 1471203v 8 - 1184301v 4 + 47385)^) c' 2 

V 9(v 4 -3) 4 (3 + v 4 ) 2 (v 8 -9) sV 

- 1458/ - 2430,,* - 2187) k - ^(3^ - 8," + 90/ + 216/ + 243)^ 2 
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